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ABSTRACT 



Aims. We investigate the stellar pancake mechanism during which a solar-type star is tidally flattened within its orbital plane passing 
close to a 10* Mq black hole. We simulate the relativistic orthogonal compression process and follow the associated shock waves 
formation. 

Methods. We consider a one-dimensional hydrodynamical stellar model moving in the relativistic gravitational field of a non-rotating 
black hole. The model is numerically solved using a Godunov-type shock-capturing source-splitting method in order to correctly 
reproduce the shock waves profiles. 

Results. Simulations confirm that the space-time curvature can induce several successive orthogonal compressions of the star which 
give rise to several strong shock waves. The shock waves finally escape from the star and repeatedly heat up the stellar surface to high 
energy values. Such a shock-heating could interestingly provide a direct observational signature of strongly disruptive star - black 
hole encounters through the emission of hard X or soft y -ray bursts. Timescales and energies of such a process are consistent with 
some observed events such as GRB 970815. 
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1. Introduction 

It has long been underlined that stars orbiting closely enough 
to massive black holes could be tidally flattened into a 
transient pawcafe-shape confi guration until full disruption 
(ICarter&Luminetl 119821 Il983h . Strongest compressions have 
been predicted to trigger thermonuclear explosion in the core of 
small main-sequence stars grazing bl ack holes from thousands 
up to millions solar masses ( Pich onj [19851: iLuminet & PichonI 
■ fT989h . The presence of specific proton-enriched chemical ele- 
' ments near galactic centres may be c ompatible with the stellar 
pancakes nucleosynthesis ( Luminet & Barbuv 1990). 

More recently, high-resolution hydrodynamical simula- 
tions have shown that the tidal compression could inde- 
pendently give rise to strong shock waves within the stel- 
[ lar matter jlC obavashi et all 120041: iBrassart & Luminet! l2008t 
iGuillochon e t al. 2008). Propagating outwards, shock waves 
may then be able to quickly heat up the stellar surface and lead 
to the emission of a new type of hard X or soft y -ray bursts. 

Aftermath of tidal disruptions of stars by massive black holes 
have already been detected from nearby otherwise quiescent 
galactic cores through giant-amplitude persistent X-UV flares 
(e.g . Bade et al. 1996; Grupe et al. 1999; Komossa & Greiner 
1999; Greiner etal. 2000; Komossa 2002; Halpern et al. 2004; 
Komossa et al. 2004; Gezariet al. 2006; Esquej et al. 2007; 
Gezari et al.1 l2008t ICappelluti et all l2009l) . These flares have 



In o ur previous article (I Brassart & Luminel"2008', hereafter 
I Paper J) , we have simulated the tidal compression orthogonal to 
the orbital plane using a one-dimensional Godunov-type shock- 
capturing hydrodynamical model, and assuming that the star 
evolved within a Newtonian external gravitational field. 

In this article, we naturally extend the model to the general 
relativistic case of a non-rotating Schwarzschild black hole and 
present some complementary results. 

In the following, relativistic equations are expressed in ge- 
ometrized units c = G - 1, greek and latin indices run respec- 
tively from to 3 and 1 to 3. 

2. Basic equations 

2.1. Orbital equations of motion and relativistic tidal field 

We consider the Schwarzschild coordinates in which the black 
hole of mass M, is described by the usual spherically symmetric 
line element 



2M.\ . I 2M. 

\dt + 1 



arised from the long term evolution when part of the liber- 
ated stellar gas was being accreted by the central black holes. 
However, the very short-timescale high-energy flares due to the 
initial compression of the star could interestingly allow to catch 
disruptions from their beginning. 



r j \ r 

+ r^ide^ + sin^edip^). (1) 

The motion of the star's centre of mass follows a timelike 
geodesic of proper time r defined by (e.g. ,Chandrasekhar.l983l) 

t^E[l-—\, (2) 
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with (■) the derivative with respect to r. Such an orbit lies in the 
equatorial plane 6 - 7r/2 and entirely depends on both constants 
of motion E and L, respectively the specific total energy and the 
specific orbital angular momentum of the star 

The tidal field applied within the star is characterized by 
the relative acceleration between infinitesimal matter elements 
viewed as a collection of test particles. Let m'' be the four- 
velocity of the star's centre of mass, and kf' the four-separation 
with a nearby test particle moving along an adjacent geodesic 



such that u^'k.j 



0. Their r elative accele ration obeys the 



geodesic deviation equation (e.g. lMisner et al 



■R 



livp 



(6) 
(7) 

with V the covariant derivative and R"^.,,, the Riemann curva- 
ture tensor Let w^^^ be an orthonormal tetrad locally defined 
at the star's centre of mass, and parallely propagated along the 
geodesic motion such that 

<o) = (8) 

^(a) 8nv = nap, (9) 

where the indice in brackets refers to a vector of the tetrad, 
and rii^v - diag(-l, 1,1,1) are the metric tensors respectively of 
the space-time around the black hole and of the flat space-time. 
For a comoving observer at the star's centre of mass, labelling 
the infinitesimal nearby matter elements by the three-position 

z' = w<;V, (10) 

the tidal acceleration ^ can thus be written in the Newtonian- 
type form (,Pirani.l956) 



(11) 



where tiie tidal tensor is defined in terms of the Riemann curva- 
ture tensor as 



Cij = -R^,p^ v/l^^ w^o) wljy 



(12) 



The symmetric traceless tidal tensor represents the first 
quadrupolar contribution of the tidal potential 



1>t = -^zVC,, + 0(|z|3). 



(13) 



The higher-order contributions relating the deviation of the star's 
orbit from a strict geodesic can be neglected considering the 
small size of the star 

It is possible to directly connect the tetrad t o the coordina tes 
of the star's centre of mass along its geodesic (lMarcklll983l) . It 
comes 



' ' "I 

vv^ij cos i/r - w^2) sin 4', 

vv^2) = "^1) sin \f/ + 'A^^2) ^'^^ 
1 
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5^ (18) 



■S'i, (19) 



with 6y the Kronecker's symbol, and where the t// variable satis- 
fies 

EL 

<A=T-72- (20) 
+ L^ 

From (fT4l)-(fT9b and the Riemann curvature tensor, the non-zero 
components of the tidal tensor (fT2l i for a Schwarzschild space- 
time take the explicit form 



Cii 
C22 
C33 



M. 



, r2_^L2 
r2 + L2 



■ sin iff 



-^(l.3^ 



r3 

1-12 T- 



1-3 ' \" C0s2|/r|, (21) 

(22) 
(23) 
(24) 

C21 = C12. (25) 

We assume that the star draws a parabolic-type orbit around 
the black hole for which E - 1 and L > AM,. When the last 
inequality arises, the black hole actually generates the gravita- 
tional barrier w hich c an reject the incoming star to infinity (e.g. 
Chandrasekha3 IT983h . From the cancellation of (O, the perias- 
tron of the orbit is located at 



r2 + L2 
-3 2 — cos ip sin ip 



AM. 



1 + 



16 



(^) 



(26) 



The critical value L - AM, dividing plunge and escape orbits 
corresponds to the minimum periastron 



' p m 



: AM, 

■ 2 r„. 



(27) 
(28) 



where = 2M, is the gravitational radius of the black hole. 
From ( l26l l, the specific orbital angular momentum defining the 
parabolic orbit is thus given by 



1/2 



(29) 



where the periastron is set through the value of the penetration 
factor 

/? = - (30) 



of the orbit within the characteristic disruption tidal radius 



M^ 



1/3 



(31) 



with Mi, and respectively the mass and radius of the star 

To solve the orbital equations ©-([Si) and ( |20] i of parabolic 
motio n with E - 1, we rewrite them as (.Frolov & Novikovl 
il998h 



1-^ 



1 1 _^J" 
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defining 



~ L 
—, 

tp - L Kg, 



1/2 

and switching to dimensionless variables 

T 

f - —, 
h 
t 

r 



(37) 

(38) 
(39) 

(40) 



(41) 
(42) 
(43) 



with (■) the derivative with respect to f. By convention f, 0, and 
(/( are set null at the periastron of the orbit, and are negative and 
positive respectively for earlier and later positions. For f > 0, 
the temporal integration of (l33T l. (l35T l. and (|36] | leads to 



. _ 2 / ^r^ -r + {f + 1- ii) 
3 \ f-^i 



+ ^ • F{^\a)- ^\+lyiE{ip\a)\ (44) 



Fiip \ a). 



iP- 1 
\h - arctan -i / r - 1 + — H ^ 



(45) 

F(^ \ a), (46) 



introducing the elliptical integrals of the first and second kind, 
respectively 



F(ip\a)= I (1 - sin^ a sin^y?)"'^^ dp, (47) 

ml/2. 



Jo 



£(<^\a)= (1 -sin^a sin^/3)"^46, (48) 



with 



(fi - arcsin 



. 2 4yU 

Sin a - 



V^^^ _ ^ + ^2 



1 +2yU 



(49) 
(50) 



The equivalent expressions for f < deduce from the previous 
ones taking into account that r is an even function of f, whereas 
(p and tfr are odd functions of f. 

2.2. Hydrodynamical equations 

Except for the definition of the tidal field, the h ydrodyn amical 
model closely follows the assumptions made in iPaper 11 . Since 
the size of the star is excessively small relative to the curva- 
ture radius of the external space-time, it is quite possible to 
keep a Newtonian description of its internal motions. The hy- 
drodynamics is described in the comoving reference frame that 



is parallely propagated along the geodesic of the star's centre 
of mass. We still restrict to simulate the one-dimensional mo- 
tion of the stellar matter along the orthogonal direction to the 
orbital plane whe n the st ar is moving within the tidal radius. 
As emphasized in Paper % such an approximation (compared to 
full three-dimensional hydrodynamics) is well justified when the 
star evolution is calculated only during the free-fall and bounce- 
expansion phases, since then the vertical motion is fully decou- 
pled from the induced motions in the orbital plane, and the latter 
are quite negligible compared to the dynamics in the vertical di- 
rection. 

Setting a vertical z axis from the centre to the poles of the 
star, the Euler's equations write 



defining 



drU + d,F(U)^S(U,z,T), 



U = {p,pv,pef , 
F{U) = [pv, P+pv^,{P+ pe)v]' 
S(l/,z,T)= [0,pg,pvg]^, 



(51) 

(52) 
(53) 
(54) 



respectively the conserved variables, the fluxes, and the sources, 
where p is the density, P the pressure, v the velocity, e the spe- 
cific total energy, and g the tidal acceleration. 
The specific total energy 



V2 

^^2 



(55) 



includes the specific internal energy s and the specific kinetic en- 
ergy. The stellar matter is considered as an ideal gas of constant 
adiabatic index 7 = 5/3 for which 



p = (r-i)pe. 

The corresponding temperature satisfies 



^ - — 



(56) 



(57) 



with the specific gas constant set for the solar composition. 
The speed of sound is locally defined by 



a = ^liy — y)yE. 



(58) 



The general three-dimensional expression of the tidal accel- 
eration is given by ( fTTT i with (I2ni-(|25]|. However, for the one- 
dimensional model, we only need to take into account the third 
component which is responsible for the compressive motion or- 
thogonal to the orbital plane. It simply writes 



g(z, t) = - 



M. 



1+3- 



(59) 



with the specific orbital angular momentum of the star defined 
by ( |29l ) and the radial coordinate deduced from the resolution of 

61. 

Finally, the conservation form of the Euler's equations ( fSTI )- 
( |54] i represents a hyperbolic system of conservation laws. It is 
solved by the Godu nov-typ e shock-capturing source-splitting al- 
gorithm detailed in lPaper iL 
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rcos(W[R.] Time[s] 

Fig. 1. Left: Parabolic-type geodesic orbit of a solar-type star 
around a 10* M© Schwarzschild black hole deeply plunging 
within the tidal radius with a penetration factor yS - 5. The solid 
circle represents the black hole's gravitational radius and the dot- 
dashed circle the tidal radius. The point on the orbit shows the 
instant of maximum compression at the centre of the star. The 
stellar orbit within the tidal radius does not significantly dif- 
fer from its Newtonian counterpart, and finally intersect outside 
the tidal radius so that the tidal field induces one compression. 
Right: Evolution of the central pressure as a function of proper 
time T. After the passage of the star through the periastron at 
T = 0, the stellar matter suddenly compresses at the centre of the 
star until the pressure reaches the magnitude needed to counter- 
act the compression's contribution of the black hole's tidal field. 
The central pressure is highest at r » 24.66. 



3. Results 

Identically to the Newtonian calculations, we are interested in 
typical encounters between a black hole of mass 
and a main-sequence star of mass M*. - Mq and radius Ri, = Rq 
taken as an initial polytrope of polytropic index 3/2. The initial 
central density, pressure, and temperature are respectively noted 
Pi,, Pi,, and Ti,. We consider cases where the star deeply plunges 
with /3 > I within the tidal radius dSTT i to simulate its resulting 
orthogonal compression. The existence of the minimum perias- 
tron ( l27b however imposes a maximum penetration factor a; 1 1 
(above this value the star enters the black hole). 

Relativistic modifications on the stellar orbit and on the 
tidal field significantly appear when the star approaches the 
black hole's gravitational radius. Due to relativistic precession, 
the parabolic-type orbit must eventually intersect once in the 
Schwarzschild space-time, and depending on the position of 
the crossing point, the star can be subjected to several suc- 
cessive compressions during its motion with in the tidal radius 
jLuminet & MarckHigSStFLaguna et alJll993h . 

Crossing orbit outside the tidal radius 

A first encounter with = 5 is illustrated on Fig.[T] The star 
is still far away from the gravitational radius for the space-time 
curvature to remain weak, and for the crossing point to be 
located outside the tidal radius. The description of the tidal 
compression process therefore proceeds as in the Newtonian 
case in three distinct phases: over-all free fall, central bounce- 
expansi on with shock waves formation, over-all expansion (see 
iPaper ll Figs. 4, 6, 7). The evolution of the hydrodynamical 
variables during the bounce-expansion phase can be followed in 

Crossing orbit inside the tidal radius 

A second encounter is illustrated on Fig. [3] for a penetration 
factor 13-9, huge enough for the star to come nearer to the black 
hole's gravitational radius and experience the space-time curva- 




Z/R. Z/R. 



Fig. 2. Density (left) and velocity (right) profiles in the positive 
vertical direction z at different proper times r during the bounce- 
expansion phase for 13 - 5. Labels stand for the following val- 
ues of T [s]: (1) 24.66, (2) 27.43, (3) 28.39, (4) 28.75. The star 
is at the periastron at t = 0. The central compression is max- 
imal at T a! 24.66. The stellar matter expands from the centre 
whereas it collapses elsewhere. Both opposite motions produce 
pressure waves that steepen into a shock wave. Powered by the 
expanding motion, the shock wave propagates outwards through 
the collapsing matter with a velocity ^ 10^ km s"', and induces a 
strong compression ratio » 3.5. The collapse stops at t ~ 28.80 
when the shock wave reaches the stellar surface a; 1 .4 s after its 
formation. 
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Fig. 3. Left: Parabolic-type geodesic orbit of a solar-type star 
around a 10* Schwarzschild black hole deeply plunging 
within the tidal radius with a penetration factor (3-9. The solid 
circle represents the black hole's gravitational radius and the dot- 
dashed circle the tidal radius. The points on the orbit show the 
instants of maximum compression at the centre of the star The 
stellar orbit winds up inside the tidal radius so that the tidal field 
induces two successive compressions. Right: Evolution of the 
central pressure as a function of proper time r. The central pres- 
sure is highest at r » -37.68 before the passage through the 
periastron, and at r ^ 89.13 after the passage through the peri- 
astron. 



ture. The stellar orbit winds up around the black hole within the 
tidal radius, which allows the tidal field to highly compress the 
star twice before and after the periastron. 

The stellar matter enters a. first supersonic free-fall phase un- 
til the central pressure highly increases to make it bounce and 
stop the compressive motion (Fig. HJi. As in the previous case, 
a strong shock wave S \ forms outwards during the bounce, and 
escapes from the stellar medium to give way to an over-all super- 
sonic expansion. However, the periastron being not yet reached, 
the tidal field becomes again prevailing relative to the central 
pressure of the star which decreases as suddenly as it previously 
increased (Fig.|5]l- 

The stellar matter is then forced to enter a second free-fall 
phase after the periastron, which occurs differently to the first 
one (Fig.|6|l. During that phase, a strong shock wave 52 forms as 
was already the case in the New tonian context for large enough 
penetration factors (see lPaper J, Figs. 8, 13). However, its forma- 
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Fig. 4. Velocity (left) and Mach number (right) profiles in the 
positive vertical dkection z at different proper times r during 
the first compression for yS = 9. The temporal evolution reads 
from top to bottom, (a): over-all free fall, (b): central bounce- 
expansion, (c): over-all expansion. Labels stand for the follow- 
ing values of t [s]: (1) -102.94, (2) -79.74, (3) -56.54, (4) - 
37.68, (5) -37.56, (6) -37.33, (7) -37.09, (8) -37.00, (9) -34.81, 
(10) -23.41, (11) -11.76, (12) 0.15, (13) 11.60. The star is at the 
periastron at r = 0. The point on the velocity profiles corre- 
sponds to the sonic point, the flow being subsonic (resp. super- 
sonic) to the left (resp. right). The homologous velocity profiles 
v(z, t) ~ z characterize the free-fall and expansion phases in the 
external gravitational field. The Mach number reaches remark- 
ably high values. The stellar surface collapses with a maximum 
velocity a; 5 x lO-' km s"', and expands with a maximum veloc- 
ity « 10'*kms The free fall stops from the centre of the star 
at T =s -37.68 when the central pressure reaches its maximum. 
The stellar matter then begins to expand while it continues to 
collapse elsewhere. Pressure waves steepen into a shock wave 
that propagates outwards with a velocity » lO^kms"' until 
it escapes from the medium at t » -36.98. It induces a strong 
compression ratio a; 3.7. The whole stellar matter then continues 
to expand. 



tion occurs quite further from the centre of the star The shock 
wave propagates slower than the stellar matter collapsing be- 
hind it, so that the latter goes through the former which is finally 
ejected from the medium without reaching the centre of the star. 
While the stellar matter continues to collapse, the central pres- 
sure begins again to increase (Fig. [T]) following the winding of 
the stellar orbit, and a third shock wave 5 3 of small amplitude 
forms. This shock wave propagates up to the centre where it col- 
lides with the symmetric shock wave propagating on the other 
side of the orbital plane. 

The reflexion of both shock waves definitely stops the cen- 
tral compression, and the stellar matter begins again to expand 
from the centre (Fig. [SJ. In front of the reflected shock wave, a 
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Fig. 5. Evolution of the tidal acceleration at the surface of the 
star as a function of proper time r. The tidal acceleration ( l59l ) 
is maximal at the stellar surface at each time. The simultaneous 
evolution of the central pressure is added for comparison. The 
star is at the periastron at r = 0. The double successive phases 
of free-fall and expansion are indicated. The decreasing of the 
tidal acceleration during the orbital motion of the star allows the 
central pressure to suddenly counteract twice the free fall, before 
and after the periastron passage. 



last shock wave 5 4 finally forms following the pressure waves 
which are generated by the opposite motions of expansion and 
collapse. This strong shock wave propagates up to the stellar 
surface and escapes followed by the previous one. The whole 
compressive motion stops to give way to a second phase of 
over-all expansion while the star moves towards the tidal radius. 

Observational signature via Xjy -ray bursts 

As already highlighted in iPaper li the shock waves gener- 
ated during the single Newtonian orthogonal compression could 
highly heat up the surface of the star in the hard X or soft y 
-ray domains. Such a mechanism may provide a direct observa- 
tional signature of strongly disruptive encounters between stars 
and massive black holes in galactic nuclei. 

The temperature profiles associated to the double relativistic 
orthogonal compression of the star are reproduced in Fig.|9l The 
different shock waves actually lead to heat up three times the 
stellar surface. A first time by the shock wave S \ formed during 
the first bounce-expansion phase. A second time by the shock 
wave 5 2 formed during the second free-fall phase, and finally 
caught up by the faster stellar surface. A third time, quite near to 
the previous one, by the shock wave S4 formed during the sec- 
ond bounce-expansion phase. At each occasion, the shock waves 
drop off a temperature a lO'^K at the stellar surface (FigfTOli. 
Therefore, it appears that the prompt emission of hard X or soft y 
-ray bursts which may be associated to the stellar pancake mech- 
anism would be interestingly intensified by the relativistic effects 
of multiple compressions. 

4. Conclusion 

We have presented first explorative results on the star's tidal 
compressions by a massive Schwarzschild black hole which take 
account of the intrinsic shock waves formation orthogonal to the 
orbital plane. Whereas the stellar core would be heated up by 
the external compressive tidal field (Fig. [TTI ). the hydrodynam- 
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Fig. 6. Density (left) and velocity (right) profiles in the pos- 
itive vertical direction z at different proper times r during the 
second free-fall phase for^S = 9. The temporal evolution reads 
from top to bottom. Labels stand for the following values of t 
[s]: (14) 56.48, (15) 68.13, (16) 72.80, (17) 76.29, (18) 78.61, 
(19) 79.31, (20) 80.69, (21) 81.40, (22) 84.31, (23) 85.49, 
(24) 87.86. The star is at the periastron at t = 0. While the 
stellar matter collapses in free fall, a shock wave 52 of strong 
compression ratio « 3.5 forms (a). The shock wave propagates 
inwards with a velocity » 10^ km s whereas the stellar surface 
collapses with a velocity 10 times superior. The former is finally 
caught up by the latter and escapes from the medium without 
reaching the centre of the star (b). The velocity profiles (c) show 
that one part of the stellar matter is slowed down going through 
the shock wave, while another part stops to collapse near the stel- 
lar surface following the ejection of the shock wave, and while 
in the front the non-affected stellar matter continues to free fall 
with a homologous velocity profile. A shock wave 5 3 of weaker 
compression ratio a; 1 .5 forms near the interface of both collaps- 
ing parts, and reaches the centre of the star atT ^ 89.13. 



ical simulations also suggest that the stellar surface would be 
heated up during the shock waves propagation. 

Besides the more restrictive scenario of the core thermonu- 
clear explosion, it is thus believed that the close passage of a star 
to a massive black hole may be directly revealed as highly ener- 
getic short-timescale bursts originating from the shocked stellar 
surface. It is well known that some observed gamma-ray light 
curves present two or more peaks with complex structures (e.g. 
lFishmanlll995l) . Could some of them be interpreted as tidally- 
induced pancake stars? Typical timescales shown in Fig. [10] for 
a star - black hole encounter along a relativistic self-crossing 
orbit are quite consistent with those of a short duration gamma- 
ray burst with two peaks of energy separated by a hundred of 
seconds, comparable to w hat has been observed e.g. in GRB 
970815 dSmith et al.ll2"002l) . The long duration GRB 060614 has 
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Fig. 7. Evolution of the central pressure as a function of proper 
time T during the second free-fall phase for j3 - 9 (enlargement 
of Fig.O. While the stellar matter suddenly compresses for the 
second time at the centre of the star, a shock wave S 3 forms at 
T a 84.31 on both sides of the orbital plane (a). Both symmetric 
shock waves propagate inwards until they collide at the centre 
of the star at t x 89.13 (b), which produces an additional (in- 
stantaneous) compression of the central matter (b) — > (c). The 
shock waves are then reflected outwards which stops the central 
compression. 




Fig. 8. Density (left) and velocity (right) profiles in the pos- 
itive vertical direction z at different proper times t during the 
second bounce-expansion phase for [3-9. The temporal evo- 
lution reads from top to bottom. Labels stand for the following 
values of t [s]: (25) 89.51, (26) 89.98, (27) 90.45, (28) 90.94, 
(29) 91.38, (30) 92.54, (31) 93.72, (32) 94.89, (33) 96.05. The 
star is at the periastron at t = 0. The point on the velocity profiles 
corresponds to the sonic point, the flow being subsonic (resp. su- 
personic) to the left (resp. right). After its reflexion at the centre 
of the star at t 89.13, the small-amplitude shock wave 5 3 
propagates outwards, in front of the stellar matter which begins 
to expand. Pressure waves induced on the still collapsing stellar 
matter leads to the formation of another shock wave 5 4 in front 
of the first one (a). This shock wave produces a strong compres- 
sion ratio X! 3.7, and reaches the stellar surface at r 96.18 (b). 
It then occurs an over-all expansion. 
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Fig. 9. Temperature profiles in the positive vertical direction z 
at different proper times t for /3 - 9. The temporal evolution 
reads from top to bottom. Labels stand for the following val- 
ues of T [s]: (4) -37.68, (5) -37.56, (6) -37.33, (7) -37.09, (8) - 
37.00, (14) 56.48, (15) 68.13, (16) 72.80, (17) 76.29, (18) 78.61, 
(19) 79.31, (20) 80.69, (21) 81.40, (22) 84.31, (23) 85.49, 
(24) 87.86, (25) 89.51, (26) 89.98, (27) 90.45, (28) 90.94, 
(29) 91.38, (30) 92.54, (31) 93.72, (32) 94.89, (33) 96.05. The 
star is at the periastron at t = 0. (a): first bounce-expansion phase 
(see Fig. |4] (b)). The shock wave S i carries out a temperature 
jump from the centre to the surface of the star and escapes from 
the medium, (b), (c), (d): second free-fall phase (see Fig.|6]l. The 
shock wave S 2 carries out a temperature jump from the surface 
to the centre of the star, but is caught up by the collapsing sur- 
face and thus heats it up before escaping. The shock wave S3 
then caiiies out a small temperature jump to the centre of the 
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Fig. 10. Evolution of the temperature as a function of proper 
time r for /3 - 9. The star is at the periastron at t = 0. The 
solid line corresponds to the temperature at the centre of the star. 
The compression induces a first maximum ^ 2 x 10^ K during 
4 s, then a second one with the same amplitude during a; 6 s. 
The dashed line corresponds to the temperature carried out along 
the vertical direction by the shock waves. The three maxima are 
reached at the stellar surface. (la): temperature increase due to 
the propagation from the centre to the surface of the shock wave 
S 1 resulting from the first bounce-expansion phase. The surface 
reaches a temperature ^ 4 x 10^ K during x 0.03 s. (Ila): tem- 
perature increase due to the propagation from the surface to the 
centre of the shock wave S 2 resulting from the second free-fall 
phase. Since it finally meets the shock wave, the stellar surface 
is heated up a: 10^ K during a 1.5 s. (lib): temperature increase 
due to the propagation from the centre to the surface of the shock 
wave S 3 resulting from the second bounce-expansion phase. The 
surface reaches a temperature a lO^K during a 3 s. It lasts 
XI 115 s between the first and second peaks of temperature at 
the surface, but a quite shorter time a 15 s between the second 
and the third ones. 



also been tentatively interpreted by o ther authors as a strongly 
disruptive star - black hole encounter dLu et al.ll2008b . 

Tidally-induced gamma-ray bursts are estimated to occur 
once every 10^-10-** years per galaxy, depending on the nuclear 
stellar density profile and the mass of the central black hole 
(IWang & Merrittll2004 . Such a rate can even be inc reased by a 
factor 10"^ in the case of a massive binary black hole (IChen et al.l 
12009 *). Since most of the galaxies - including our own Milky 
Way - harbour a central massive black hole, and since the full 
observable universe is transparent in X or 7 -ray domains, several 
events of this kind would be then detectable each year 

To well characterize this signature and make valuable com- 
parisons with observational data of high-energy bursts, simula- 
tions have to consider the full stellar pancake's structure and 
the radiative transfer of energy through the layers of the star. 
From this perspective, three-dimensional high-resolution shock- 
capturing simulations have been recently performed in the con- 
text of a Newtonian gravitational field, and have shown that 
shock waves actually form in several other directions within the 
star with the consequence o f isotropizing the surface tempera- 
ture (iGuillochon et al.ll2008l) . 
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Fig. 11. Evolutions of the maximum central density (top), max- 
imum central temperature (middle), and duration during which 
the highest values are maintained (bottom) as functions of the 
penetration factor /? in the case y - 5/3. The rightside up (resp. 
upside down) triangles indicate the first (resp. second) compres- 
sion. The points show for comparison the single compression 
calculated for a Newtonian gravitational field, and the dashed 
lines correspond to the power laws predicted by the affine model 
(TCarter & Luminetl 119831: see also IPaper L Eqs. (39)-(41) and 
Fig. 16). In bottom figure, the characteristic internal timescale 
of the Starr* = (1/Gp*)'/^ 



